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Abstract 

A variety of unitary gauges for perturbation theory in a background field is considered 
in order to find those most suitable for a Hamiltonian treatment of the system. We select 
two convenient gauges and derive the propagators for gluonic quantum fluctuations 
immersed in background configurations. The first one is a unitary generalization of the 
usual Coulomb gauge in QED which preserves the decoupling of two propagating polar- 
izations from the instantaneous one. The second possibility is the axial light cone gauge 
which remains ghost free also in the presence of a background. Applications of the for- 
malism to the spectrum and dynamics of QCD at the confinement scale, such as hybrid 
states, are briefly discussed. 
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1 Introduction 



The identification of relevant degrees of freedom fT], |, H in the non-perturbative region of QCD 
is the crucial step towards a consistent theory of hadrons. One of the ways to disentangle 
this problem in an economical way is to separate in Euclidean space first non-trivial 'vacuum' 
configurations B^, e.g. to split the total gluonic field into a background and a quantum 
fluctuation 



A^ = B^ + (1.1) 

where the field is primarily responsible for the long range gluonic correlation functions. 
The qualitative picture for S^-fields is the ensemble of lumps described by their collective 
coordinates. For the self consistency of the separation (p..l|), these coflective coordinates should 
not be distorted after inclusion of quantum fiuctuations a^. 

Recently, this idea has been emphasized |], with the suggestion to consider perturbation 
theory in a confining background which leads to the area law asymptotics of the averaged 
Wilson loop. The (path) integral over the background fields can be reformulated as a 
summation over irreducible correlators of vacuum field strength tensors by means of the cluster 
expansion. In this way, the existence of fiux tubes between quarks is related to decaying 
correlators with particular Lorentz structures 0. Although the procedure appears reasonable, 
and is supported by the results from the simpler 2 + 1 compact U{1) theory P| where the 
monopole background fields generate at large distances the 'frozen' string, no consistent theory 
exists for such background fields in QCD. At present, one must postulate that the dynamics 
of the vacuum configurations leads to a finite string tension, which can be expressed as the 
infinite sum over contributions from the vacuum correlators of all orders. 

Since the confining background cures all infrared singularities [Q, it enables one to investi- 
gate the dynamics of the a^-perturbations at large distances as well. The important point is 
that the a^-fiuctuations immersed into such vacuum fields describe the effective long distance 
excitations of the string 0, |^ frozen at the level of the vacuum contribution. They should 
be considered, besides quarks, as the relevant degrees of freedom at the scale of confinement. 
For instance, the interpretation of the one a^-gluon exchange amplitude (averaged over the 
5^-fields) of Fig. la requires introducing the first excitation of the string as in Fig. lb. In 
general, the QCD string between quarks will be divided by valence gluons into a corresponding 
number of elementary 'frozen' pieces. 

Within this picture, the basic question is how to describe the gluonic excitations ('con- 
stituent gluons') interacting via the background in a nonperturbative way. Because of gauge 
invariance, in eq. (|l.l|) contains (unphysical) pure gauge components. Also the presence of 
non-decoupled propagating ghosts inherent in the standard formulation |^, 0, Q of perturbation 
theory in a background field with the usual gauge condition 



D^{B)a^' = (1.2) 

will obscure the analysis of bound states. It is the main goal of this paper to explore these 
problems and to select unitary background gauges where ghosts are either absent or nonprop- 
agating. This will allow us to describe bound states in terms of physical polarizations of 
only. 

To apply the propagators D^^{B) =< a^a^ > to a given problem, we suggest to use the 
multichannel Hamiltonian approach P, |^ recently generalized [|lOl to this kind of fluctuating 
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strings. It provides a scheme to derive from the propagator both the diagonal elements of the 
Hamiltonian for states with a fixed number of propagating constituent gluons as well as the 
nondiagonal terms which mix different Fock states. In this way, we can investigate glueballs or 
hybrids. We stress that Dn^{B) provides us with both spin-dependent and spin- independent 
nonperturbative interactions between a^- gluons and other constituents at the scale of con- 
finement. We note that the diagonal approximation for such states based on the standard 
nonunitary background gauge (|L^) (with the propagating ghosts) has recently been suggested 
in ref. 0. 

The plan of the paper is as follows. In the next section, we discuss further the special role 
of the unitary background gauges for the analysis of QCD bound states. Then, in section 3, 
we review the derivations of D^^, in the abelian case within the path integral approach. 

In section 4, the standard analysis of perturbation theory in background fields is extended 
to prove the well known background gauge invariance 0, |l], H for the case of the special unitary 
gauges to be described in sections 5 and 6. There, the explicit form of the propagators D^i,{B) 
will be derived. Finally section 7 contains a brief sketch of the applications and the conclusions. 



2 The special role of unitary background gauges 

Before constructing explicitly the propagators in the unitary background gauges, we discuss in 
more detail why they are important for the analysis of gluonic excitations in bound states. 

As stressed above, the key point is that if dynamical gluons appear as constituents of a 
bound state, such as in glueballs or in hybrid states, non-physical polarizations of D^^y ( which 
occur in a nonunitary gauge) obscure both the interpretation of the quantum numbers of a 
hadron and the dynamical scheme of interactions between the valence constituents. Due to 
gauge invariance, the contributions from 'pure gauge' polarizations should drop out (together 
with ghost contributions) in the final result. But in order to invoke physical intuition, it is 
important to have the transparent and economical formulation which is provided by unitary 
gauges where ghosts are either absent or nonpropagating. Then all independent propagating 
polarizations in D^^, can be interpreted as physical degrees of freedom. 

This special role is illustrated already at the level of standard perturbation theory if a bound 
state is involved. In the calculation of deep inelastic scattering amplitudes, it is the axial or 



planar gauge (see |TT[ for the refs.) in which only the planar Feynman graphs contribute 
to the leading logarithmic asymptotics which allows for a selfconsistent introduction of the 
hadronic structure functions. In such a gauge, there is a direct correspondence between properly 
chosen independent fields and the "physical" degrees of freedom. This conceptual necessity of 
a particular gauge stands in contrast to the usual perturbative analysis where the expansion 
is restricted to a given order in the coupling constant which makes results explicitly gauge 
invariant. In that case, the choice of gauge is essentially a technical matter. 

Another example of a dynamical scheme where gluons become 'constituents' is the ap- 
proaches based on the light-cone QCD-Hamiltonian |p. As in any quantum field theory, a 
multichannel Hamiltonian must be introduced in order to reproduce, via iterations, the ampli- 
tudes which initially are expanded in time ordered Feynman graphs. In the absence of ghosts, 
wave functions of Fock states with extra dynamical gluons will describe only the "least" number 
of physical polarizations. In this approach the unitary light cone gauge (for the standard 
perturbation theory) insures this property. 

To apply effectively the multichannel Hamiltonian approach to the dynamics of gluonic ex- 
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citations in a background we will therefore look for gauges where 

i) ghosts are either absent or do not propagate, 

ii) there is an unambiguous and simple way to separate in D^y{B) the non-propagating polar- 
ization from the propagating ones. 

A gauge satisfying condition i) will be called unitary; it reproduces correctly the number 
of independent physical propagating polarizations interacting with the background fields. The 
second requirement gives an additional selection of unitary gauges leading, as we will see, to 
the most economical use of the Hamiltonian approach. 

Gauge invariance obviously allows to eliminate one of the four components of the a^-gluons 
in a background. But since they are generally not 'on mass-shell', it is by no means evident how 
to formulate a generalized transversality condition by which only two propagating polarizations 
are retained and which is equivalent to the non-propagation of ghosts. We will show in this 
work that there are two most suitable unitary gauges and derive the corresponding propagators 

The first one generalizes the ordinary Coulomb gauge propagator of QED and is based on 
the condition 

N,{B)a, = {D,{B) + 2D^\B)Foi{B))a, = 0. (2.1) 
leading to a propagator D^^{B) which satisfies 

Do.(S) = 0, z = 1,2,3. (2.2) 

Here D^^{B) and -F^^ = /"''^F^^(5) are the standard covariant derivative and field strength 
tensor respectively. 

The Dqq{B) part corresponds to the instantaneous Coulomb exchange modified by the 
interaction with the background. As for the spatial components Dik{B), we will show that the 
generalized projector 

p-b = s^,5''' - {N,{B)Nr'iB)NkiB)r'' (2.3) 

allows one to represent Dik{B) in the sandwiched form 

D,k{B) = PUB){K{By^)mnPnk{B). (2.4) 

This gives the correct transversality condition for the propagating physical components orthog- 
onal to the instantaneous polarization in accordance with eq. ( p.2| ). 

All other unitary modifications of the Coulomb gauge in the presence of a background do 
not admit the decoupling condition ( |2.2D important for the further applications of D^^. As a 
result, requirement ii) is not satisfied. 

The second unitary gauge can be directly obtained from the light cone gauge constraint 

{n^a^Y = 0, nl = 0. (2.5) 

Other axial gauges (with 7^ 0) being unitary in the presence of background fields do not 
satisfy condition ii). 

The transversality condition of eq. ( |2.4| ) for the nonperturbatively interacting gluons has 
two immediate consequences. First it gives selection rules for the allowed quantum numbers 
of states with a fixed number of gluonic constituents. A well known example is the (Landau- 
Yang) theorem for the forbidden total momentum of the system of two real photons. Similarly, 
this condition decreases the dimensions of multiplets with fixed allowed quantum numbers. 
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3 Path integral derivation of QED propagator 
in Coulomb and axial gauges 

Let us first review the procedure of gauge fixing performed directly in tlie patli integral repre- 
sentation for tfie generating functional Z{J). 

We will formulate the criterion for a gauge to be ghost free in a way which can be easily 
generalized to the perturbation theory in a nonabelian background. Our starting point is the 
standard expression in Euclidean space 



Z{J) = J DA^lc.F.exp - J d'x (^^ 



and we impose the gauge invariance condition 



df^J^ = 0. 



To evaluate Z{J) in the form 
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Z[J) ~ exp 

one must separate explicitly in eq.( p.l|) the gauge zero modes A^f^ 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



of the quadratic term or, in other words, fix a gauge. The standard way is to insert into 
eq. (|3.1|) the Faddeev-Popov unity in the form 



1 = det 



W5 

6u 



(3.5) 



where A^^ = + d^g and f{A^) is the gauge fixing function satisfying f{A^) = 0. Thus 

1 



Z{J) 



*(/(AJ)exp 



d X 



4^7 



2 ^lJ.v "I" Jf^A^ 



(3.6) 



where we omitted the volume of gauge modes / Dg. 

By fixing the gauge in this way, the number of independent polarizations in D^,^ is al- 
ways decreased from four to three. There remains the freedom to select one non-propagating 
polarization. 

Rewriting F^^ via the kinetic matrix K^i, as 

and taking for simplicity zero sources = 0, one gets 



Z(0) 



det^ 



5/ nl/2 



_detKik_ 

Here, i, k refer to the three components of Kf^^, selected by 6{f{A)). 



(3.7) 
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Instead of counting independent propagating polarizations of D^jy, the absence of ghosts 
can conveniently be formulated as follows: detKn^ should have a factorized form 




detKii, = c^et I — j det {PimK^nP, 
d Sf 



d{dl) 6u 



(3.8) 



where the operator (|5)^^ must be non-propagating and the corresponding determinant in 
eq.(|3.7[) must cancel exactly the ghost factor in the numerator of eq. (|377|). We use the 
notation det in order to stress that the last determinant in eq. ( p.8|) denotes formally the result 
of the integration over the two propagating physical polarizations (selected by the projector 
Pik {PimPmk = Pik)) ■ They determine completely (if eqs. (|3.8|) are satisfied) the statistical sum 
Z(0) 

Z{0) = \(M (PimKmnPnk)]'^'^ ■ (3.9) 



In QED, the ghosts decouple always and the first condition of eq. (|3.8| ) holds; only the 
second one is nontrivial. In the presence of a nonabelian background -B^, the ghosts interact 
with the i?^-fields and factorization (|3.8|) of the determinant in the statistical sum Z{0, B) is not 
insured for every gauge. This is related to the fact that in gaussian approximation Z^'^\Q, B) 
of the statistical sum is not invariant under the choice of background gauges for quantum 
fluctuations. But a background gauge is nevertheless unitary (ghosts do not propagate) if the 
second condition of eq. ( p.8| ) is met. 

In the remainder of this section our aim is to gain experience in how to use the freedom in 
gauge fixing for bringing the propagator D^^{x — y) into the form the most convenient for the 
Hamiltonian technique. 

A) Coulomb gauge 

Instead of imposing from the begining the gauge condition ( |3.5| ) in the form 



f{A,) = d,A, (3.10) 
one can start from the 'physical' conditions and require that D^j, satisfies 

d 



9^92)^00 id,) = 0, D^of{d,) = D'r^\d,) = 0, z = 1,2,3 (3.11) 

which we would like to preserve in the presence of a background. The second condition of 
eqs. (|3.11|) implies that the spacial and temporal components, Ai and Aq, are decoupled, while 
the first one insures that Aq does not propagate in time (and is responsible for the instantaneous 
interaction) . 

Since there are only two transverse polarizations for ordinary photons, one can impose an 
additional constraint on Dik{d^) via insertion of the gauge unity ( p.5| ) in a form which leads 
to eqs.( p.llD . To understand which function f{A^) brings D^^, into the required form, we first 



write j d'^x {jF^^ as 

j d^x^F^, = ^J d^xi-Aod^Ao + 2AMAi - {AdlS.kAk - AAdkAj,}). (3.12) 
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Comparing with ( |3.11| ), we see that the mixed term AododiAi should be canceled. This imme- 



diately leads to the standard form (|3.1CI|) of f{A^ 

Summarizing, we see that starting from conditions one can selfconsistently determine 

the form of gauge fixing function f{A) leading to the ghost free gauge. In Section 4 we will 
generalize this idea to the case of perturbation theory in a nonabelian background. 

The main steps leading to Z{J) which will be repeated in the nonabelian case are as follows. 
In order to get rid of 6{diAi), we introduce the projector on transverse states 



didk 
■'i 

Decomposing into transverse and longitudinal parts 



P^k=\5ik-^y = 0. (3.13) 



A = PikAk + (1 - P)ikAk = Aj- + Al (3.14) 

one obtains 

6id,A,) = (rfet(-92))-i/2<5((l - P),kAk). (3.15) 
Therefore, the integration over (1 — P)ikAk = A| is eliminated and Z{J) takes the form 



Z(J) = j DAoDAi{det{-d'^)y/^exp +\ j d^x ("4(^5^0+ 



(3.16) 



Because only the transverse components enter, the resulting propagator will have the sand- 
wiched form 

Dik = PilDlmPmk, (3-17) 

given after integration over Aq and Aj- by the final representation ( p.3| ) for Z{J) 



ZiJ) = {deti-dl)r'exp 



^ d^x (jiPii-g^PikJk + Jq-^Jq 



(3.18) 



We stress that the remnant of the ghost determinant in eq.( p.l6| ) has been cancelled by the 
integration over the instantaneous component ^40. The expression for the statistical sum (3.9) 



Z(0) = (rfet(-aj))^i 

is to be interpreted as that for two (transverse) propagating polarizations. 



(3.19) 



B) Axial gauge with ^ 

In the abelian case (or in the case of standard nonabelian perturbation expansion) there is 
another conventional ghost free gauge in which the propagator includes only two propagat- 
ing polarizations, the well known axial gauge [nA) = 0. The gauge fixing function can be 
conveniently choosen as 



f{A) = {nd){nA). 



(3.20) 
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so that the second condition in eq. ( p.8|) , written in terms of hght cone variables, holds. Here, 
n is an arbitrary 4-vector. If it is light-like, the corresponding gauge is the light cone gauge. 

All are conceptually equivalent to the elementary cases n = (1,0) and n = (0, 1, 0_l). 
The derivation of the propagator is similar in both; we discuss only 



n^ = (l,0). (3.21) 

(3.22) 

^k. (3.23) 

and the projector Pik is defined by eq. ( p.l3[ ). The representation ( |3.23| ) for Kik provides 
the separation of nonpropagating and propagating polarizations such that integration over the 
former cancels the ghost factor det^/^(— 5q). This separation is analogous to that of eq. ( ^.16] ); 
represents the nonpropagating component. 

Rewriting DAi as DA^DA", we arrive at the final representation 



In full analogy with the previous subsection, one obtains 

Z{J) = j DAi{dei{-dl)f'^e^v +^ J d^xi^AKi^Af, - 2J,A,) 

where 



Z(J) = (det(-a2))-iexp 



3_ 
2 



J d^x U{P^P)^,J, + J,((l - P)l(l - P)),,J, 



. (3.24) 



In accordance with gauge invariance, Z{0) has the same form (and interpretation) as in the 
Coulomb gauge. We note that in D^^, the nonpropagating part represented by the last term in 
the exponent of eq. ( p.24| ) replaces the instantaneous Coulomb exchange part of eq. ( |3.18| ). 

C) Light cone gauge, = 

To complete the analysis, we consider also the light cone axial gauge [nA) = 



n^ = (3.25) 

where the light cone coordinates for a 4-vector are b± = -^{bz ±bo) and we choose the gauge 
fixing function / as before. 

This choice of n requires a Minkowski metric and straightforward calculations lead to the 
following analogue of eqs.( |3.16D and ( p.22| ) 



Z{J) = I DA^DA^{det{-dl)y/^exp J d^x (^-^A^K^.A, + 2J,A?j (3.26) 



where fx, v run over +, — , ±= {i, A;} and 



(3.27) 



The form ( |3.27| ) of K^^ implies that the integration over the non-propagating component A_ 
exactly cancels the ghost factor ((iet(— (9i))^/^ with the result 



Z(J) = J DA^exp d'x{[AK,kAk + 2g^Mi] + 



+ {g\Jn) - {A^d^){nd)}—-^{g'{Jn) - {nd){d^A^)} 

[noy 



(3.28) 



Consequently, the final form of Z{J) reads 



Z{J) = {deti-dl))-'exp 



+ — / d X 



{Jn) 



{ndy 



iJn) - 



which corresponds to 



D 



-9f^u + 



(nd) 



dl + ie ' 



(3.29) 



(3.30) 



The sum in the exponent of eq. (|3.29| ) naturally provides us with the separation of propagator 
(3.30) into the instantaneous 

(3.31) 



D 



inst 



and propagating parts 



r)Prop 



9^u + 



(nd) 



{ndy 



dl + ie' 



(3.32) 



This will be useful when we extend our considerations to the nonabelian background. It is easy 
to check that the transversality condition is satisfied 



iV^D^r = 0, N, = d^^+m,{nd) 



with the new 4-vector m„ 



(3.33) 
(3.34) 



1, m_ = m_L = 

which implies that there are only two propagating polarizations as before. 

Summarizing, we see that eliminating in Z{J) one polarization by inserting the gauge unity 
(|3.5| ), one generates in all three cases a corresponding ghost determinant which is exactly 
canceled by the integration over the retained physical 'nonpropagating' polarization. In other 
words, the determinant 

det{K^,) (3.35) 

of the kinetic matrix for the three retained polarizations assumes a factorized form ( p.8[ ) where 
the part which cancels the ghost determinant can be separated from the one describing two 
propagating physical polarizations. 



4 Gauge fixing procedure for the perturbation expan- 
sion in a nonabelian background 



As in the abelian case, we begin with the path integral form |^, |T^| of the generating functional 
in the presence of a background 



Z{J,B) 



Da^ det 



Sf{B^,a'^ 
5uj 



5(/(5,a))exp 



d X 



^{F^^{B + a)y + JXu 



MM 



4.11 
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where f{B^, '^'^)Lo ~ keeping in mind unitary gauges. In what follows the standard notations 



12| are used: 

F^,{B + a) = F^.T'' = F^,{B) + D^{B)a, - D,{B)a^ - z[a^, a,] (4.2) 



with 



Df = 8^5"^ + r^^B^ = 9^5"* - ary^Bl = df,6^^ - i[B^,y^ (4.3) 



[O'fii (^u\ 



alal[T\r]=ir^T^alal (4.4) 
For zero sources = 0, the action density ~ F'^^{B + a) obeys the gauge symmetry 

{B^ + a^) ^ U{B^ + a^ + id^)U-^ (4.5) 

which in the infinitesimal form reads 

B^ + a^^ B^ + a^ + D^{B + a)uj, (4.6) 
and where D^{A)uj can be expressed according to eq.(4.3) as 

T'^id^.u^ + r^'u^A') = T^d^uj" + (cj X AY). (4.7) 



There are two related issues which distinguish the gauge fixing procedure in eq.( f4.1| ) from 
the one for the standard perturbation theory with no background. 

The first one is the well known fact []^ that in the presence of B^ it is possible to select for 
the a^-fields a gauge which preserves the invariance of Z{J, B) under gauge transformations of 



Bfj^ and the sources 



B^-^U{B^ + id^)U-\ J^^UJ^U-\ (4.8) 
Usually this invariance is exploited within the background field method J^, [1^ to simplify 



the renormalization scheme and multiloop calculations. Here we would like to stress the inte- 
gration over background configurations resulting in particular properties of the asymptotics 
for averaged Wilson loops. To build up the formalism in terms of gauge invariant objects like 
Wilson loops (with possible spin- insertions, see for example [0), it is necessary to represent 
physical amplitudes in terms of gauge invariant combinations of Greens functions in the back- 
ground fields such as < ^f(a;)\l/(y) > and < a^(x)a^(?/) > for the the matter and the a^-fields, 
respectively. Consequently, it is economical to work with such a gauge for in which the 
Greens function D^y{x,y\B) =< a^{x)ay{ii) > is transformed homogeneously with respect to 
background gauge transformation 



D^,{x,y\B) ^ U{x)D,,{x,y\B)U~\y). (4.9) 



This property follows if Z{J, B) is invariant under 
The second aspect is the freedom to split the gauge variations of B^ + field in eight ways 



between B^ and fields: 



S^a^ = D^{B + a)uj, 5'^^'^B^ = (4.10) 
(5(2)a^ = D^{B)u, 6^^^B^ = {u x a^) (4.11) 
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= (u;xa^), b^^^B^ = 


D^{B)u 


(4.12) 




= D^{a)u, d^^^B^ = 


{uj X B^) 


(4.13) 




= {u;xB^), 6^'^B^ = 


-- D^{a)u 


(4.14) 




d^uj, 6^^^B^ = {ux 


(B^ + a^)) 


(4.15) 



5(^)a^ = {ujx {B^ + a^)), 6'^'^B^ = d^u (4.16) 

where we omit the eighth irrelevant sphtting S^^^a^ = 0. 

Since no sufficiently comprehensive discussion of these points exists in the literature [T^ , 
we will present the detailed analysis required to construct the nonstandard unitary gauges in 
the next two sections. 

First, we single out the splittings leading to a Z{J, B) which is invariant under ( ^.8| ). Since 
the action of eq. (|4.1j ) itself is invariant under all seven forms of the variation (the third one, 
(4.12), corresponding to that of eq. ( |4.8|) ), possible non invariances must come from the gauge 
fixing determinant det ^ — -^^^^ '° ^ of the Faddeev-Popov unity in eq. (|4.1| ). Thus, in order 

to insure invariance under ([4.8| ), must transform homogeneously under the infinitesimal 
variation (|4.12|) (with infinitesimal parameter a)), i.e. 



(4.17) 



For simplicity, we ffist take covariant gauges and consider the function f{B, a) (linear in a^) 
written in terms tranforming homogeneously under ( 4.12|) such as 



D^{B)a,, F^,{B)a„.... (4.18) 

This implies that in all interesting cases i?^ enters only inside covariant derivatives D^{B) or 
field strength tensors F^^{B). 

Consequently, the partial derivatives ^-^^^'"^ and ^^^f^'""^ also consist of structures which are 
transformed homogeneously. As a result, if gauge variations 

6uj 6B Suj 5a 6u 

are to satisfy eq. ( ^4.17] ), the quantities and must be also constructed from combina- 
tions transforming themselves in accordance with this equation 0. 

With this in mind, we condclude that only the three ffist variations (^4.10|) - ([4.12|) lead to 



an invariant generating functional Z{J,B). In all other cases, and consist of 'wrong' 



(in the above sense) combinations of 5^-fields. 



2 



It is not difficult to insure in all interesting cases that indeed there is no "fine tuning" between the two 



terms of eq.(4.19) and each should satisfy the condition separately. 
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To illustrate this point, we work out the standard background gauge case 

r{B,a) = {D,{B)a,r = 



Combining with eqs.( |4.1U| ) 



5r_ 

6a^ 



ab 



D^{B + a 



5^'^B, = 



(4.20) 
(4.21) 

(4.22) 



one gets the standard result || 



In a similar way, eqs.( [4.11| ) leading to 



{D,{B)D''{B + a)) 



ab 



(4.23) 



yield 



ab 



{D^{B)D,,{B) + a^a^) 



ab 



while eqs.( [4.12D 



c„c\ab 



6u^ 



implies the final expression 



ab 



~t[D,{B),a, 



D,{B 



lab 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



Simple but tedious calculations confirm our general conclusion that these three forms of gauge 
fixing of the a^- field do exhaust all possibilites to maintain the invariance of Z{J, B). 

We observe that the determinants of the operators ( [4.23[ ) and ( [4. 251 ) in the leading order 



0) are reduced to 



det{-Dl{B)), 



(4.28) 



which generalizes the abelian factor ( p.l9[ ) and represents two propagating ghost polarizations 
interacting with the background. Note that the determinant of the function ( |4.27| ) vanishes 
formally for = 0. In the following, we will always use the first kind of gauge splitting ( [4.10|) . 
In Landau gauge, eq. ( [4.20| ), one obtains in leading order 



Z(2)(0,S) 



det{-Dl{B)) 



11/2 



'det{-P,p{B){5,^Dl{B) - 2Fp,{B))PUB)) 



(4.29) 



where P^y = 5^^ — D^^D^ and the determinant in the denominator indicates the formal 
integration over three components of selected by condition ([4.20|) (see the next section for a 
similar detailed derivation). We immediately conclude that the conditions of eq. (p.8D are not 
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satisfied. Ghosts are propagating and in D^j, all three retained polarizations are propagating. 
Obviously, this gauge is not unitary. 

From the above analysis it is clear that the same arguments can be applied to noncovari- 
ant gauges to maintain the invariance under ( [4.81) . We must only insure that gauge fixing 
function f{B, a) is constructed from elements with entering via the (noncovariant) combi- 
nations of long derivatives and field strength tensors. This is exactly what we will exploit when 
constructing the nonabelian background generalizations of Coulomb and axial gauges. 



5 Coulomb background gauge 

We first apply the above considerations to obtain the most suitable unitary generalization of 
the Coulomb gauge for the perturbation theory in a background field which satisfies eq. ( p . 1 1 ) . 
We will also derive the corresponding propagator D'^J;{B) which is required for applications. 

As in the abelian case, we start with the quadratic approximation for F^^{B + a) to write 
the generating functional Z^'^^ in the form 



Z^^\J,B) = J Da^det 
where 



^W/(^^,a^) 



UifiB, a)) exp +1 J d'x (^-^X^K;t{B)al - 2JX^ 



= [D%B)5,, - D,{B)D^{B) - 2F,^{B)X 



and 



f: 



ab _ rabc rpc 



We neglected in the expansion of F^j^{B + a) the contribution of the term linear in 

4<D-(i?)F;,(S) 



(5.1) 
(5.2) 

(5.3) 

L 

(5.4) 



assuming that it merely renormalizes the parameters of the effective action F'^^{B) expressed 
through the proper collective coordinates. This is realized, for instance, in the 1 + 1 gas of 
kinks and antikinks ||T^. We stress that this assumption is necessary for the existence of a 
selfconsistent separation (|1.1|) of the field into the two parts. 

In order to obtain the gauge fixing function /(.B, a) leading to conditions ( p.ll| ) we must 
cancel the mixed terms in a^K^ya^, 



aoKoiai = -ao£'o(A + ^Foijai 



which leads to the following generalization of the abelian condition (|3.1CI|) 



riB, a) = {D,{B) + 2Do i(S)Fo,(S))^V = {N,{B) 



(5.5) 



(5.6) 



The operator Nf'{B) selects the two transverse 'physical' polarizations. This gauge fixing 
function yields the second requirement (p.8|) . 

We note that the gauge fixing function f{B, a) of eq. ( |5.6| ) can not be obtained by a simple 
substitution of di by Di{B) in the abelian form (|3.1CI|) . 
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With this f{B,a), it is straighforward to insert the corresponding gauge unity into the 
generating functional; we obtain 



Z^'^\J,B) = J Da^det{-Ni{B)Di{B + a))5{Ni{B)ai) : 



X exp 



where 



Kf^iB) = {6MB) - 2F,k{B) - D,{B)D,{B)r . 



In order to ehminate 6{Niai) we introduce, as in the abehan case, the projector 
P,kiB) = (^d,k - N,{B)j;^^Nk{B)^ , N,P,k = PM = 0, 



(5.7) 
(5.8) 

(5.9) 



singhng out two physical propogating polarizations among the three components. With 

ai = Pikak + (1 - P)ikak = aj- + aj (5.10) 

and 



S{Nia,) = {deti~Nl))-'/^5ia!), 



(5.11) 



one arrives at 



r.r2W T I ^deti-NiiB)DAB + a)) 



deH\-Ni{B)) 
+ ai{K,k + NiNk)ai} - 2Joao - 2J,^ai 



(5.12) 



where we substituted Kik Kik = Kik + NiNk for a later convenience (note that Niaf = 0). 

Also here, only the transverse current = PikJk enters and thus Dik will have the sand- 
wiched form 



Afc(S) = P,m{B){k{B)-').^r.Pnk{B), 



(5.13) 



generalizing the abelian expression ( p.l7|) . Consequently after integration over ao,aj- one ob- 
tains the final representation for Z^'^\j, B) 



Z^^\J,B) 



det{-Ni{B)Di{B + a)) 
[det{-NKB))det{-Dfm] 



1/2 



{dct{^ PimKmnPnk)) 



-1/2 



X 



X exp 



Q f 

-— J {JiDikJk + JoDooJo)d^x 



with 



Doo{x,y\B) =< x\ 



Dm 



\y >, Doi = Dio = 



and 



D,k{x,y\B) =< x\Pij{B){6jiDl - 2Fji - D,Di + N,Ni)-^ Pik{B)\y > 
The two equations (|5.15|) ,( prT^ ) are the central result of our paper. 



(5.14) 

(5.15) 
(5.16) 



14 



Let us now give a short interpretation of these expressions. We conclude from eq.( |5.14]) that 



ghosts are present but do not propagate in time (the factor in curly brakets of eq. ( ^.14|) is equal 
to unity when = 0). Indeed the ghost kinetic term (NiDi) differs from the corresponding 
one (Df) for the nonpropagating aQ-component only by the spin dependent part (the last 
term in eq.(p.6|)). In contrast to eqs.( ^4.28| ),( |4.29| ) ghost are not propagating and enter only 
the loop corrections to the instantaneous Coulomb exchange given by Dqq at the tree level. 
Consequently, this is a unitary gauge and ghosts do not introduce independent propagating 
degrees of freedom (described by Dik) which is sufficient for the desired economical formulation 
of the Hamiltonian approach. 

We point out that the expression ( |5.16| ) for Di^ gives the correct form for both spin- 



independent and spin-dependent interactions of transverse gluons (decoupled due to conditions 
(|5.15|) from the instantaneous component) with the background. 

If the conditions ( p. 11 ) are not imposed, there is a variety of unitary gauges generalizing 



the abelian Coulomb one ( |3.1CI| ). Indeed every gauge where Ni{B) has the same limit 

N.{B)-^d, for B^^O (5.17) 

like Ni{B) and obeys the second condition of eq. ( ^.8] ) will lead to nonpropagating ghosts. 
Because of the coupling between the and the two retained Oj-components (-Doi 7^ 0) the 
choice of the nonpropagating component is not ambigous in this case. It is the simple form 
(|5.15|) and ( |5.16| ) of D^i, which renders the choice ( p.6|) of Ni[B) the preferable one for further 



applications with the multichannel Hamiltonian. 



6 Axial background gauges 

Next, we derive Z{J, B) in quadratic approximation when the gauge condition is imposed in 
the form 

{n^a^r = 0. (6.18) 

As before, the term linear in (see eq. ( |5.4|) ) is neglected also here. 

First we analyse whether ghosts are also absent in these gauges if there is a background. 
For this purpose it is sufficient to consider Z*^^)(0, B) and look for a representation in the form 
of eq.( p.9| ). It will be shown that only the light cone (n^ = 0) axial gauge remains ghost free in 
the presence of an arbitrary background. Again we will obtain the propagator and give a short 
interpretation. 



A. The case ^ 

As in the abelian case, for 7^ it is sufficient to consider n = (1,0) or n = (0, 1, 0_l). Both 
elementary choices are similar, we take only the temporal gauge 

n=(l,0), (6.19) 

for which in Euclidean space 

Z^^\0,B) = I Daidet"\-Dl{B))ei^^[+^j d^xaiK,u{B)ak], (6.20) 
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with 

= [(DliB) + DfiB))6,, - D,{B)D,{B) - 2F,^{B)Y\ (6.21) 

We note that in this gauge the ghost factor (iet^/^(— Dq(5)) does not depend on the a^-fields. 
Equation ( |6.20D immediately leads to 



Z(2)(0,S) 



det{-Dl{B)y 
det{-Kik{B)) 



(6.22) 



Again we see that Z^'^\{),B) is different for different background gauges because the gaussian 
approximation to Z(0, B) implicit in perturbation theory is not invariant under a change of the 
background gauges. This raises the question whether there is a 'best' gauge. 

As in the abelian case, the absence of ghosts is equivalent to the existence of a projector 

P^^{B) = 5'^'5.,-lyMmj^^M,{B)\i , P' = P 



= M„{Mkak) + Pikak (6.23) 



1 

such that the following generalization of the first condition in ( p.8| ) holds (the second one is 
obviously satisfied) 

det{-Kik) = det{-Dl)d^t{-P,^Kl^Pnk). (6.24) 

In other words, ghosts are absent if there exists an operator M^^{B) such that the result of 
integration over the instantaneous component defined as 6 = [M^ak) [Dai = DbD{Pa)j_) 



cancels exactly the ghost factor. The last condition results in 

MiDlSikMk = MiKikMk (6.25) 
which means in turn that the 3-dimentional kinetic matrix Kn. = — ^ikD^ 

= DfSik - DiDk - 2F,k (6.26) 

has a continuum of zero modes I Mi }-— h ) . This finally leads to the constraint 



det[-kik] = (6.27) 

reproducing the condition ( |6.25| ). 

To demonstrate that eq. (|6.25|) can not be satisfied without an extra constraint on F^^{B) 
let us consider first the simpler 2 + 1 case where 

det[-kik] ~ det{-Dl + {DiD2 + F2i)Df{D2D, + Fu)). (6.28) 

The operator in eq. ( p.28| ) can be represented as 

(FuD^'Di - DiD^'Fu) + FuD^^Fu, (6.29) 

and one concludes that it vanishes in general only if 

[Di,F,k]=0, i,k = 1,2. (6.30) 
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This is not surprising because Kik of eqs.( |6.26|) , ( |6.28D coincides with the kinetic operator (|5.2|) 
for the 1 + 1 case and therefore has gauge zero modes DiOj [jl], in the presence of a background 
which satisfies the classical 1 + 1 equations of motion ( |6.30 ). In the 2 + 1 case at hand, any 
classical background satisfies 

[A,^^fc] = -[I^o,^ofc] (6.31) 
which is equivalent to eq. ( |6.30| ) only if in addition 



(6.32) 



is imposed. 

It is not difficult to show that in 3 + 1 dimensions, eqs. (|6.3CI|) with i,k = 1,2,3 are necessary 
to satisfy eq. (|6.25|) . Therefore, for 7^ it is impossible to cancel the ghost factor exactly for 
an arbitrary background field, even if it is only classical. Thus we conclude that in general the 
background axial gauge ( |6.18|) is not ghost-free when 7^ 0. 

Still, for all gauges where Mf{B) satisfies the second condition (Op and 



Mm 



M,{B) 



--di, Bu 







one gets 



Z^^\0,B) 



' det{-Dl)det{-Ml 
det{MiKikMk 



1/2 



det '^'[{PK'P),k] 



(6.33) 



(6.34) 



and the factor in curly brakets becomes unity for = when we recover the abelian form 
( |3.19| ). Consequently, the only propag ating part of Z^^\0,B) arises from the integration over 
{Pikdk) and ghosts lead merely to loop corrections to the nonpropagating contribution of the 
exchange. As a result, all such gauges ( |6.33| ) can be called unitary. 

The impossibility to cancel det^^'^{—DQ) corresponds to the ambiguity in the definition of the 
projector Pik onto transverse propagating polarizations since all Mi{B) of eq.( |6.33| ) are suitable 
(if the second condition ( |3.8| ) is not violated). Consequently, one can not select unambigously 
the non-propagating component from the rest. 



B. Light cone gauge = 

In this gauge with 

n_ = 1, = nf = 0, 
the quadratic approximation to Z{J, B) is given by 



Z^^\J,B) = J Da_Da,{det{-Dl{B))Y'^ x 

+^ / d'^^^-^o.^lK^,u{B)a^ + 2(J+a_ + Jiai)| 



where 



xexp 



K++ = -Dl, = -D+D, - 2F+„ = -D,D+ - 2F,+ , 
Kik = Dl6ik — DiDk — 2Fik 



(6.35) 



(6.36) 



(6.37) 



and k stand for the two perpendicular components. 
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From these equations one obtains after integration over the a_ component 



X exp 



Dai exp 



2^2 



j d'^x{aiKikak - 2g^ai{Ki+K^\jj^ - Jj)} 



X 



(6.38) 



where 

kt^ = Ktl: - {K,^{K^^)-'K^ur. (6.39) 

We see that the ghost factor in eq. (6.36) is cancelled exactly by the integration over a_ as in 
the abelian case. The absence of ghosts in the light cone background gauge makes it the most 
convenient one (among the axial gauges) for the Hamiltonian approach built up in terms of 
physical polarizations alone. 

Integrating out the a^-fields, one obtains Z'^'^\j, B) in the representation of eq. (|3.3D 



Z(2)(J,S) = {detA',fc)}-iexp|-^ j d^x{L,{J)k-^^ Lj,{J) + J+{K++)-\j^ 



(6.40) 



where 



L,(J) = J,,-K,+ (K++)-V+ (6.41) 

The physical interpretation of eq. ( |6.40|) is similar to the abelian case (|3.29|) . The determi- 
nant in the preexponent reproduces the statistical sum Z(0, B) in term of two propagating 
polarizations. In the exponent itself, the first term gives the part of the propagator D^J'°^[B) 
responsible for these two propagating polarizations. The second term supphes us with the 
modified instantaneous exchange 



D^-:\B) 



DUB) <"^' 

generalizing the abelian expression ( p.31| ).We note that the choice of physical propagating 
polarizations is unambiguous. 

Similarly to Dik{B) of eq.( p.l(: ) in the Coulomb background gauge, the piece D^J°^{B) given 

by 

Drpr^{B) = -{g^^ - n,{nKnr\nK)f)K;^{g^y - {Kn)^{nKn)-\y) (6.43) 

provides, in particular, the light cone representation for the spin-dependent interactions of 
physical polarizations with the background. We note that it is more complicated than in the 
Coulomb case. 



7 Applications of unitary gauges 

As pointed out in the introduction, the basic motivation of our work is to design a dynamical 
scheme to investigate QCD bound states including gluonic excitations at the scale of confine- 
ment ("constituent gluons"). In this section we sketch how the propagators can be used to this 
aim. 

In order to simplify the discussion and to separate the confining dynamics we are interested 
in here from the effects of chiral symmetry breaking for light quarks, we will consider the case of 



spinless quarks only. For heavy quarks the spin- dependent interactions can be related ||T3[ via 
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a cluster expansion to the same irreducuble correlators of background fields which constitute 
the averaged Wilson loop. For light quarks these interactions are connected to chiral symmetry 
breaking and an additional information on the confining configurations is required. Work in 
this direction is now in progress. 

The Feynman-Schwinger representation allows one to write the Greens function of the qq 
state in a path integral form. For spinless quarks in the quenched approximation we have 

G = j dsiDzids2Dz2exp[-Ki - K2] < W{C) >B+a (7.1) 
where < W{C) >b^^ is the Wilson loop operator averaged over all gluonic fields = + 

< WiC) >B+a= \ I DBDadet f 



Af{B,a))x 



X exp 



J F'^j^{B + a)d^x Pexp i j {B^ + afj,)dx 



(7.2) 



and the integration over DB^ implies a summation over all relevant collective coordinates ex- 
cluded from Dttfj^ in the standard way (we omit for simplicity the corresponding orthogonality 
conditions ensured via the Faddeev-Popov trick). The contour C consists of quark and anti- 
quark trajectories zi, Z2 and Ki and K2 are the standard quark kinetic terms whose explict 
form will not be required here. 

Quantum fluctuations can be taken into account by expanding the path ordered exponent 
Pexp[i JciBf^i + a^)dXf^] in powers of a^. The result can be represented as 

00 „ 

< W{C) >B+a=< W{C) >B + E ^" / {W^'"\C; x{l), x{n))) dx{l) . . . dx{n) (7.3) 

n=l 

where < W^^'^ >b is a Wilson loop with n insertions of the fields [Q. For example, the 
quantity 



j < W^'^^ >B dx{l)dx{2) 



= 9 
+Oig') 



dx,il)dx,{2){^"c^^ixil),x{2)\B)D'^Z'^^i^^^^^ 



B 



+ 



(7.4) 



(with D^^'^^ = t^"t^^D^|j,) describes (apart from the 0{g^) corrections) the one gluon exchange 
in Fig. la. Here $(7. denotes necessary parallel transporters in the fundamental representation 
along the subpathes Ci of the initial contour C 



<^>cAx,y\B) = Pexp 



{BT)dx, 



(7.5) 



while D^iy is the propagator of the fluctuations in the background. 

We emphasize that the subscript B in eqs. ( [7.3|) , ( [7.4|) refers to retaining only the 5^-field 
in all relevant path ordered exponents. The averaging procedure is still performed with the 
action density F^^{B + a) as in eq. ( |7.2| ). 

To incorporate < W^^^ (C) > b consistently into the bound state formalism, a multichannel 
Hamiltonian approach has been proposed [|lOl . Within this approach, time ordered diagrams as 
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in eq.( [7.3|) can be reproduced by iterations of the effective Hamiltonian. The problem is thus 
reformulated as that of a relativistic many body system with mixed channels. 

Under general circumstances, a selfconsistent application of the multichannel Hamiltonian 
formalism to a quantum field theory requires the light cone frame (coordinates). For this 
purpose, the covariant generalization |T5| of the Coulomb gauge ( |3.10| ) can be used where the 



new propagator takes the form (|5.15|) , ( ^.16|) in the co- moving Minkowski coordinates, obtained 



by a Lorentz transformation in accordance with the total velocity of the system 

<^ = ^M = 7(l,^,0x), n« = 7(^,1, 0±), n(f'^) = (0,0,nx). (7.6) 



One can express n^^^ via the conventional light cone coordinates to recover the conceptually 
important suppression of pair creation from the vacuum. In the limit V ^ 1, the new Dqq 
part expressed in terms of light cone coordinates corresponds to the instantaneous Coulomb 
exchange boosted from the rest frame to the infinite momentum one |T^. Consequently, it 
makes sense to work with D'^°\B) in the form of eqs.( |5.l5D , ( |5.16| ), which is, in addition. 



somewhat simpler than D^y{B) of eqs.( |6.42D , (|6.43| ) in the light cone gauge. Keeping this in 
mind, we are led to the following strategy |]TU|. 

Take bound states with light quarks where the long distance dynamics clearly dominates. 
Then one can treat in Euclidean space the hard and soft parts of the a^-field separately as 
fast and slow subsystems. First one integrates the hard part with four-momentum squared 
> ^ where Tg is the scale at which the area law of < W{C) >b starts and the effective 
string between gluons and quarks forms. 

Apart from the inducing new effective vertices, this also leads to a renormalization of the 
running coupling constant at the scale ^ where it freezes due to the confining vacuum 

With this effective action the soft part (with Euclidean four-momentum squared < t^-) of 
the a^-field is handled, together with the 5^-field, in the following way. 

At the initial step one must take into account the contributions from the first term in 
eq.( |7.3D , < W{C) >b, plus the sum over all orders in the instantaneous exchange DQQ{x,y\B) ~ 
(5(x+ — y+) from the rest. In the following, we will always replace the integration over the 
background fields by the sum over irreducible correlators of the field strength F^u{B) with 
the help of the cluster expansion . At the hadronic scale, this enables one to use the expansion 
in where Tg ~ 0.2 — 0.4 fm is the decay length of the correlators and < r > is the 
caracteristic distance between constituents. In this way one retains ^ the first few local terms 
(such as the area and perimeter law contributions) for quantities like the Wilson loop averaged 
over 5^-fields only 

<W{C)>B=exp{-aS + pi + ...). (7.7) 

where S,l,a and p are the area, length, string tension and effective mass respectively. We note 
that at the level of the vacuum average, the contribution of the higher derivative structures 
(such as the curvature) is relatively suppressed for large smooth countours. 

Denoting this part of < W{C) >B+a as < W^'^^'^C) >B+a, one can introduce the effective 
qq Lagrangian defined formally by 

L{qq)dt+ = - In j dsiDzids2Dz2 exp[-is:i - K2] < W^^'^\C) >B+a (7.8) 

After continuing to Minkowski space, the Legendre transformation gives us the Hamiltonian 
H{qq) for the two valence quark Fock sector \l/(g^). 
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We emphasize that L{qq) contains on the hght cone both the frozen flux tube contribution 
171 arising from the area law asymptotics (|7.7|) of < W{C) >b and also the soft part of the 
instantaneous Coulomb exchange determined by Dqq{B) which describes the interactions with 
the vacuum fields. Diagramatically, it corresponds to Fig. 2. 

At the next stage, one can calculate those time ordered diagrams of < W{C) >B+a which 
are obtained by the iterations of the the two-coupled channel Hamiltonian 



H 



acting on a Fock state vector 



(2) 



H{qq) V{qq qqg) 

V{qqg qq) H{qqg) 



^>{qq) 
^{qqg) 



which includes also the state with one valence gluon (see Fig. lb). For the case 
approximation to H has been considered in ref. 



(7.9) 



(7.10) 



this 
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Generally, i/*^^-* must be evaluated according to the scheme which we discuss at the example 
of the simplest subset of the Dik iterations in Fig. 3. Time intervals with the propagating 
gluons are described by the Hamiltonian H{qqg) of the first excitation of the string 'frozen' 
during the rest of the time. It can be determined already from the simple diagram of Fig. 4a 
(where the valence gluon propagates during the entire time) in the same way as H{qq). 
This diagonal ansatz, for an arbitrary number of gluons has been suggested recently in 
We note that in order to calculate H{qqg) fully, one must take into account all instantaneous 
exchanges between valence quarks and the gluon (as in Fig. 4b). 

Without these corrections, the elements of H{qqg) are given by the Legendre transformation 
of the Lagrangian matrix Lik{qqg) (continued to Minkowski space) which is formally given by 



/ Lik{qqg)dt+ = - In [dsiDzids2Dz2 exp[-Ki - K2] < <l>(x, y\B)Dik{x, y\B)^{y, x\B) >b] 
Jo 

(7.11) 

corresponding to Fig. 4a. Apart from the kinetic and the spin-dependent terms, there are at 
large distances two frozen strings (Fig. lb) connecting the a^-gluon with the quark and anti- 
quark, respectively. They appear explicitely after using the Feynman-Schwinger representation 
of Dik{B) which implies a path ordered exponent in the adjoint representation along the gluon 
trajectories over which the path integral is performed 0. Cluster expansion techniques 
allow to reformulate the averaging in eq. ( [7.11 ) as an average for two Wilson loops in the 
fundamental representation with the closed contours (Ci + Cg), (Ci — Cg) bounded by the cor- 
responding quark and gluon trajectories. As a result, the area law asymptotics for these Wilson 
loops < W{Ci ± Cg) >B induces at large distances a (frozen) string between the gluon and 
each quark. The sandwiched form (eq. ( |5.13| )) insures that there are only two transverse prop- 
agating polarizations. Note that eq. (|7.11|) relates Dik of eq. ( ^.16| ) to both spin-independent 
and spin-dependent effective interactions of the valence gluon. 

Non-diagonal elements V{qqg —>■ qq) can be obtained uniquely from the requirement that 
the amplitude of Fig. la is reproduced in the second iteration of H^"^^ (see ref. |l^ for details). 

For the case B^ = 0, the Hamiltonian H^'^^ would lead [0 to the existence of a continuum 
part of the spectrum. Here the inclusion of the confining background allows to invoke (on light 
cone in particular) the confining QCD string and its excitations which is still not elaborated 



in the standard Hamiltonian approach p (see |]19l for a discussion). 
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In principle, one can continue this procedure to include step by step higher order Fock states 
which are relevant for a problem at hands. There are some indications JlOl that there exists, at 



least for the low lying bound states a dynamical parameter (in addition to the moderate value of 
coupling constant a^t frozen due to the confining background) which justifies this expansion. 
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Figure captions 

Fig. la One a^-gluon exchange contribution to < W{C) >B+a- 

Fig. lb Modes of the QCD string arising (on hght cone) in the amphtude of Fig. la. 

Fig.2 The part of < W{C) >B+a responsible for the dynamics in the valence quark Fock 
sector. 

Fig. 3 The simplest subset of exchanges from the part of < W{C) >B+a relevant for the two 
Fock states approximation (7.9) to the Hamiltonian. 

Fig. 4a The simplest amplitude for reconstruction of H{qqg). 

Fig.4b Amplitudes with the instantaneous exchanges leading to a corrected H{qqg). 
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